Method of Determining the Yaw Rate of a Target Vehicle

ABSTRACT

This disclosure describes a radar system configured to estimate a yaw-rate and an over-the-ground (OTG) velocity of extended targets in real-time based on raw radar detections. This disclosure further describes techniques for determining instantaneous values of lateral velocity, longitudinal velocity, and yaw rate of points of a rigid body in a radar field-of-view (FOV) of the radar system.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 16/004,487, filed Jun. 11, 2018, which claims priority to European Patent Application No. 17175548.1, filed Jun. 12, 2017, the disclosures of which are hereby incorporated by reference in their entireties herein.

TECHNICAL FIELD OF INVENTION

This disclosure generally relates to a radar system suitable for an automated vehicle, and more particularly relates to a system that calculates a yaw-rate of the target. This invention relates to a vehicle on-board method of estimation of planar motion parameters of objects (targets such as another vehicle) detected by a (e.g. Doppler) radar. It has particular application to detected the planar motion parameters of a vehicle by a host vehicle equipped with a radar system.

BACKGROUND OF INVENTION

Aspects of the invention relate to an improved method of determining the instantaneous values of lateral velocity, longitudinal velocity and yaw rate of any point of a rigid body in the radar field-of-view (FOV). Aspects are applicable (but not limited) to the estimation of vehicles for automotive perception systems and can be used in Active Safety, Driver Assistance and Autonomous Driving applications.

When an automated or autonomous host-vehicle is preceded by a target-vehicle traveling forward of the host-vehicle, it is advantageous for the system that controls the operation (e.g. steering, brakes, engine) of the host-vehicle to have knowledge of the yaw-rate of the target vehicle. Knowledge of the yaw-rate of the target vehicle can be useful to, for example, temporarily allow a closer following distance because the target-vehicle is turning out of the travel path of the host-vehicle.

In prior art, no algorithm was reported for an instantaneous estimation of full planar motion of rigid body target objects based on raw detections of a single radar.

SUMMARY OF THE INVENTION

In one aspect is provided a method of determining the yaw rate (ŵ_(t)) of a target vehicle in a horizontal plane by a host vehicle equipped with a radar system, said radar system including a radar sensor unit adapted to receive signals emitted from said host vehicle by said target, comprising: emitting a radar signal at a single time-point instance and determining from a plurality (m) of point radar detections measurements captured from said target vehicle by said radar sensor unit in said single radar measurement instance, the values for each point detection of range, azimuth and range rate; [r_(i), θ_(i), {dot over (r)}_(i)]; determining the values of the longitudinal and lateral components of the range rate equation of the target (c_(t), s_(t)) from the results ({dot over (r)}_(i), θ_(i),) of step a); where the range rate equation is

${{\overset{.}{r}}_{i,{cmp}} = {\left\lbrack {\cos\theta_{i}\ \sin\theta_{i}} \right\rbrack\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}},$

where {dot over (r)}_(i,cmp)={dot over (r)}_(i)+u_(s) Cos θ_(i)+v_(s) sin θ_(i) where u_(s) is the sensor unit or host vehicle longitudinal velocity and v_(s) is the sensor unit or host vehicle lateral velocity; determining the orientation angle of the target (γ_(t,scs)); d) determining the target center (x_(t) and y_(t)) from the results (r_(i), θ_(i)) of step a); e) determining a line l_(PH) perpendicular to the orientation of the target and passing through the center of the target (x_(t,c,scs), y_(t,c,scs)) from the results of step (c) and (d); f) determining a line l_(CA) passing through the center of rotation of said target and the position of radar sensor unit of said vehicle; from the steps of b); g) determining the intersection point of the lines l_(CA) and l_(PH) from steps e) and f) being the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] of the target; and h) estimating and yaw rate {circumflex over (ω)}_(t) from the position of the center of rotation found in step g) and components of range rate equation of the target (c_(t) or s_(t)) of step b).

Step c) may comprise determining the orientation angle of the target (γ_(t,scs)); from the values of range and azimuth (r_(i), θ_(i),) of said point detections.

Step c) may comprise determining the orientation angle of the target (γ_(t,scs)); from L fit, Hough transform or rotating caliper methodology from the results of step a).

Step b) may include estimating components of range rate equation of the target (ct, st) said point detection measurements of azimuth and range rate [θ_(i), {dot over (r)}_(i)].

Said estimating may comprise applying a Least Square method.

Said line l_(PH) of step e) may be determined and defined from the following equation:

$y = {{{{- \frac{1}{\tan\gamma_{t,{scs}}}}x} + {\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + y_{t,c,{scs}}} = {{a_{1}x} + {a_{0}.}}}$

Said line l_(CA) of step f) may be determined and defined from the following equation:

$y = {{{- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}x} = {{b_{1}x} + {b_{0}.}}}$

In step g) the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] of the target may be determined from the following equations:

${{\overset{\hat{}}{x}}_{t,{COR},{scs}} = {\frac{b_{0} - a_{0}}{a_{1} - b_{1}}{and}}},{{\overset{\hat{}}{y}}_{t,{COR},{scs}} = {\frac{{a_{1}b_{0}} - {a_{0}b_{1}}}{a_{1} - b_{1}}.}}$

In step h) the yaw rate may be determined from the following equations:

${\hat{\omega}}_{t} = {{\frac{{\overset{\sim}{s}}_{t,{scs}}}{- {\overset{\hat{}}{x}}_{t,{COR},{scs}}}{or}{\overset{\hat{}}{\omega}}_{c}} = {\frac{{\overset{\sim}{c}}_{t,{scs}}}{{\overset{\hat{}}{y}}_{t,{COR},{scs}}}.}}$

The method may include additionally determining estimates of longitudinal velocity û_(t,i), lateral velocity {circumflex over (v)}_(t,i) of certain target point from the value of yaw rate and the coordinates of the center of rotation of the target (y_(t,COR,wcs)x_(t,COR,wcs)) using the following equation:

$\begin{bmatrix} {\hat{u}}_{t,i,{scs}} \\ {\hat{v}}_{t,i,{scs}} \end{bmatrix} = {\begin{bmatrix} {\left( {{\overset{\hat{}}{y}}_{t,{COR},{scs}} - y_{t,i,{scs}}} \right){\overset{\hat{}}{\omega}}_{t}} \\ {\left( {x_{t,i,{scs}} - {\overset{\hat{}}{x}}_{t,{COR},{scs}}} \right){\overset{\hat{}}{\omega}}_{t}} \end{bmatrix}.}$

BRIEF DESCRIPTION OF DRAWINGS

The present invention is now described by way of example with reference to the accompanying drawings in which:

FIG. 1 shows a target co-ordinate system;

FIG. 2 shows a vehicle coordinate system;

FIG. 3 shows a sensor coordinate system;

FIG. 4 illustrates how to calculate velocity vectors at the locations of three raw detections;

FIG. 5 which shows the result of the cloud algorithm for yawing targets;

FIG. 6 shows an example of the determination of target orientation based on spatial distribution of raw detections;

FIG. 7 illustrates a geometrical interpretation of one method;

FIG. 8 shows an example of test results of a maneuver;

FIGS. 9 a to d illustrate the movement relating to FIG. 8 ;

FIG. 10 shows an example of test results of a maneuver;

FIG. 11 a to d illustrate the movement relating to FIG. 10 ;

FIG. 12shows an example of test results of a maneuver; and

FIGS. 13 a to d illustrate the movement relating to FIG. 12 .

DETAILED DESCRIPTION

Reference will now be made in detail to embodiments, examples of which are illustrated in the accompanying drawings. In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the various described embodiments. However, it will be apparent to one of ordinary skill in the art that the various described embodiments may be practiced without these specific details. In other instances, well-known methods, procedures, components, circuits, and networks have not been described in detail so as not to unnecessarily obscure aspects of the embodiments.

‘One or more’ includes a function being performed by one element, a function being performed by more than one element, e.g., in a distributed fashion, several functions being performed by one element, several functions being performed by several elements, or any combination of the above.

It will also be understood that, although the terms first, second, etc. are, in some instances, used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first contact could be termed a second contact, and, similarly, a second contact could be termed a first contact, without departing from the scope of the various described embodiments. The first contact and the second contact are both contacts, but they are not the same contact.

The terminology used in the description of the various described embodiments herein is for describing embodiments only and is not intended to be limiting. As used in the description of the various described embodiments and the appended claims, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will also be understood that the term “and/or” as used herein refers to and encompasses all possible combinations of one or more of the associated listed items. It will be further understood that the terms “includes,” “including,” “comprises,” and/or “comprising,” when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.

As used herein, the term “if” is, optionally, construed to mean “when” or “upon” or “in response to determining” or “in response to detecting,” depending on the context. Similarly, the phrase “if it is determined” or “if [a stated condition or event] is detected” is, optionally, construed to mean “upon determining” or “in response to determining” or “upon detecting [the stated condition or event]” or “in response to detecting [the stated condition or event],” depending on the context.

Accurate estimation of the yaw-rate and over-the-ground (OTG) velocity is important for many driving-assistance systems. Described herein is a radar system configured to estimate the yaw-rate and OTG velocity of extended targets (largely, for vehicle tracking) in real-time based on raw radar detections (i.e., range, range-rate, and azimuth). As used herein, the term ‘extended-targets’ is used to refer to targets that present multiple, spaced-apart scattering-points so the term ‘extended-target’is understood to mean that the target has some physical size. The various scattering-points are not necessarily individually tracked from one radar scan to the next, so the number of scatter-points can be a different quantity and/or each scattering point have a different location on the extended-target in successive radar scans.

The invention determines instantaneous values of lateral velocity, longitudinal velocity and yaw rate of any point of a rigid body (such as another vehicle) in the radar field-of-view (FOV). Generally, a host vehicle is equipped with a radar system where reflected radar signals (detection) from another vehicle in the field of view are processed to provide data in order to ascertain these parameters. In order to do this various conditions and requirements are needed. The target (rigid body/vehicle) needs to be a distributed target, i.e. provide a plurality of detections from the same target.

Also assumed is an approximation of the distributed target by a rigid body model which is e.g. appropriate for vehicles (passenger cars, trucks, motorbikes, trains, trams, etc.), though not generally applicable to vulnerable road users.

Radar detections received by the host vehicle form the target provide raw data provide data with respect to the position of the radar transmit/receive element/unit e.g. the Cartesian position of the detections or the Polar coordinates (azimuth angle, range). By using e.g. Doppler techniques, the range rate can also be determined.

In the subsequent concept description, the following conventions and definitions are used:

World Coordinate System

As is convention an inertial coordinate system with the origin fixed to a point in space is used—it is assumed the co-ordinate system does not move and does not rotate. Conventionally the coordinate system is right-handed; the Y-axis orthogonal to the X-axis, pointing to the right; the Z-axis pointing into the page and positive rotation is to the right of the X-axis; see FIG. 1 which shows such a co-ordinate system with origin 1 and a non-ego vehicle (target) 2.

Vehicle Coordinate System

The origin may be located at the center of the front bumper 3 of the host vehicle 4 as shown by FIG. 2 . The X-axis is parallel to the longitudinal axis of the vehicle. The coordinate system is right-handed with the Y-axis orthogonal to the X-axis, pointing to the right, the Z-axis pointing into the page and positive rotation to the right of the X-axis

Sensor Coordinate System

Origin located at the center of the sensor unit/radome. The X-axis is perpendicular to the sensor radome, pointing away from the radome. The coordinate system is right-handed: Y-axis orthogonal to the X-axis, pointing to the right; Z-axis pointing into the page; Positive rotation to the right of the X-axis. FIG. 3 shows a sensor origin 5.

In aspects of the invention and with prior art techniques, the velocity and the yaw rate of the host vehicle is assumed known. The host over the ground (OTG) velocity vector is defined as:

V _(h)=[u _(h) v _(h)]^(T),

where u_(h)—host longitudinal velocity and v_(h)—host lateral velocity.

Sensor mounting position and boresight angle in the vehicle coordinate system are also assumed known; the following notations are used: x_(s,VCS)—sensor mounting position, longitudinal coordinate; y_(s,VCS)—sensor mounting position, lateral coordinate; and γ_(s,VCS)—sensor boresight angle.

The sensor(s) Over the Ground (OTG) velocities are assumed known (determined from host vehicle motion and sensor mounting positions). Sensor velocity vector is defined as V_(s)=[u_(s) v_(s)]^(T) with u_(s)—sensor longitudinal velocity and v_(s)—sensor lateral velocity. At each radar measurement instance, the radar unit/sensor captures m raw detections from the target. Each raw detection is described by the following parameters expressed in the sensor coordinate system: r_(i)—range (or radial distance), θ_(i)—azimuth angle, and {dot over (r)}_(i)—raw range rate (or radial velocity) i=1, . . . , m.

Target planar motion is described by the Target over-the-ground velocity vector at the location of each raw detection: V_(t,i)=[u_(t,i) v_(t,i)]^(T), where: u_(t,i)—longitudinal velocity at the location of i-th raw detection, and v_(t,i)—lateral velocity at the location of i-th raw detection.

Target planar motion can be described as well by: V_(t,COR)=[ω_(t) x_(t,COR) y_(t,COR)]^(T) where ω_(t)—target yaw rate, x_(t,COR)—longitudinal coordinate of the center of target's rotation, y_(t,COR)—lateral coordinate of the center of target's rotation, x_(t,c)—longitudinal coordinate of the center of target's bounding box, y_(t,c)—lateral coordinate of the center of target's bounding box, and γ_(t)—orientation of the target object

FIG. 4 illustrates how to calculate velocity vectors at the locations of three raw detections (depicted by reference numeral 6) captured from the same rigid body target and the yaw rate of that target. Center of the target's rotation depicted with reference numeral 7. The range rate equation for a single raw detection is given as follows: {dot over (r)}_(i)+u_(s) cos θ_(i)+v_(s) sin θ_(i)=u_(t,i) cos θ_(i)+v_(t,i) sin θ_(i).

To simplify the notation, the notion of a compensated range rate is introduced and defined as: {dot over (r)}_(i,cmp)={dot over (r)}_(i)+u_(s) cos θ_(i)+v_(s) sin θ_(i), with {dot over (r)}_(i,cmp)—range rate compensated of i-th raw detection.

Then the equation is reduced to: {dot over (r)}_(i,cmp)=u_(t,i) cos θ_(i)+v_(t,i) sin θ_(i).

Range rate equation in vector form

${\overset{˙}{r}}_{i,{cmp}} = {{\left\lbrack {\cos\theta_{i}\ \sin\theta_{i}} \right\rbrack\begin{bmatrix} u_{t,i} \\ v_{t,i} \end{bmatrix}}.}$

Range rate equation in general vector form

${{\overset{˙}{r}}_{i,{cmp}} = {\left\lbrack {\cos\theta_{i}\ \sin\theta_{i}} \right\rbrack\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}},$

with: c_(t)—longitudinal range rate equation solution, and s_(t)—lateral range rate equation solution.

Velocity profile is used as range rate equation solution synonym. Note estimated values are denoted with a hat and Least Square solutions are denoted with a tilde.

Problem Formulation

The problem to be solved can be phrased as follows: estimate velocity vectors {circumflex over (V)}_(t,i) and yaw rate {circumflex over (ω)}_(t) using raw detection measurements [r_(i), θ_(i), {dot over (r)}_(i)] captured from a rigid body target in a single radar measurement instance. Since the locations of the three raw detections are known (by means of direct measurements and sensor mounting position), the equivalent problem formulation is: estimate the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] and yaw rate {circumflex over (ω)}_(t) using raw detection measurements [r_(i), θ_(i), {dot over (r)}_(i)] captured from a rigid body target in a single radar measurement instance.

Cloud Algorithm

Previously the case of a straight-line moving distributed target was considered. This restriction simplifies the estimation problem as the velocity vectors at the location of each raw detections are identical, i.e.: V_(t,i)=[u_(t,i) v_(t,i)]^(T)=[u_(t) v_(t)]^(T)=V_(t), for i=1, . . . , m.

The Cloud Algorithm (CA) was proposed to estimate over-the-ground lateral v_(t) and longitudinal u_(t) velocity of the “cloud” of detections coming from the same target. This was achieved by Least Square solution to the problem defined as follows:

${\overset{.}{r}}_{i,} = {{\left\lbrack {\cos\theta_{i}\ \sin\theta_{i}} \right\rbrack\begin{bmatrix} u_{t} \\ v_{t} \end{bmatrix}}.}$

The algorithm proved to be a very effective technique for instantaneous estimation of target velocity. Additionally, under the additional assumption of a negligible sideslip, the angle {circumflex over (γ)}_(t)=tan⁻¹({tilde over (v)}_(t)/ũ_(t)) can be used as an estimate of the targets heading.

In D. Kellner, M. Barjenbruch, K. Dietmayer, J. Klappstein, and J. Dickmann, “Instantaneous lateral velocity estimation of a vehicle using Doppler radar,” in Proceedings of 16th International Conference on Information Fusion, Istanbul, Turkey, 2013. The same problem and the same theoretical basis for the estimation of lateral velocity of a straight line moving object was considered. Enhancement to the Cloud Algorithm was made by means of executing RANSAC algorithm to identify outliers and executing orthogonal distance regression (ODR) to solve error-in-variables problem for the modified formulation of the original problem. This approach improved robustness of the solution in comparison to the original Cloud Algorithm solution. Computational complexity and the requirement to solve an optimization problem are the major drawbacks of the proposed approach, especially when an application in a production embedded system is to be considered.

Cloud Algorithm Solution for Yawing Targets

Previously the application of the cloud algorithm to the estimation of target's motion without the restriction on straight-line path was investigated. Such situation is shown in FIG. 4 . The over-the-ground velocity vectors at the location of each detection are determined as follows:

$\begin{bmatrix} u_{t,i,{wcs}} \\ v_{t,i,{wcs}} \end{bmatrix} = {\begin{bmatrix} {\left( {{\overset{\hat{}}{y}}_{t,{COR},{wcs}} - y_{t,i,{wcs}}} \right)\omega_{t}} \\ {\left( {x_{t,i,{wcs}} - {\overset{\hat{}}{x}}_{t,{COR},{wcs}}} \right)\omega_{t}} \end{bmatrix}.}$

The range rate equation for each raw detection was derived to be: {dot over (r)}_(i,cmp)=(t_(t,COR,wcs)−y_(s,wcs))ω_(t) cos θ_(i)+(x_(s,wcs)−x_(t,COR,wcs))ω_(t) sin θ_(i).

It was then shown that the Least Square solution to this problem results in: ũ_(t)=(y_(t,COR,wcs)−y_(s,wcs))ω_(t), and {tilde over (v)}_(t)=(x_(s,wcs)−x_(t,coR,wcs))ω_(t). One therefore achieves: ũ_(t)(x_(s,wcs)−x_(t,COR,wcs))={tilde over (v)}_(t)(y_(t,coR,wcs)−y_(s,wcs))

This is a result of a major importance. One can conclude (conclusion 1) that by considering range rates and azimuths from a cloud of detections coming from the same rigid body target and captured in a single look of a single radar, it is possible to estimate the position of the center of rotation and the yaw rate. Further (conclusion 2) it is possible to determine the line at which the center of rotation lies. This line passes through the sensor. This is shown in FIG. 5 which shows the result of the Cloud Algorithm for yawing targets. Calculations based on the three raw detections denoted 6 results in the line (of the center of detection) 8 passing through the sensor and the center of rotation 7. The exact location of the center of rotation 7 not specified by Cloud Algorithms.

In the results highlighted above are followed by the proposition of estimating the full motion of the target using detections captured in multiple looks (instances) of a single radar. Such an estimation is possible under the assumption of a constant yaw rate during the time of these multiple looks. It is also necessary for the host vehicle to be moving. Geometrical interpretation of this method is that one needs to find a (stationary) point of intersection of two lines passing through the (moving) sensor in multiple time instances. Because angles of the two lines are likely to be very similar (the faster the host moves the greater the angle), the approach is prone to inaccuracies.

In D. Kellner, M. Barjenbruch, J. Klappstein, jurgen Dickmann, and K. Dietmayer, “Instantaneous full-motion estiamtion of arbitrary objects using dual doppler radar,” in Proceedings of Intelligent Vehicles Symposium (IV), Dearborn, Mich., USA, 2014, the same problem was considered and the solution to the Conclusion 1 was to take measurements from two sensors. The Conclusion 2 was not mentioned by the authors, but one supposes e that the lines passing through both sensors would intersect in the center of target's rotation. The authors then applied RANSAC and ODR techniques previously presented in order to improve the robustness of the solution.

The drawbacks of the solution proposed above include that the target needs to be in the field-of-view of multiple sensors (in automotive radar configurations, common parts of field-of-view of multiple radars are usually small). For good accuracy and robustness of the solution a significant distance between the two sensors is desired. The scope of such installation is significantly limited by the dimensions of the host vehicle. One either needs to assume constant yaw rate between the looks of two radars or to synchronize the measurement instance of both sensors. Additional data processing applied to make the solution more robust are characterized by significant computational complexity.

Yaw Cloud Algorithm

The concepts reviewed so far have not assumed any motion model of the target. An Extended Kalman Filter has been formulated based on a constant turn motion model and measurement model derived from a cloud algorithm. The algorithm is called YCA (Yaw Cloud Algorithm). In the proposed formulation measurements from a single radar are considered, but the possibility of extension to multiple sensors is also highlighted.

One object in the current invention in comparison to the above, is not to use time-filtering, and not assume any underlying motion model for the target. The approach does not suffer from difficulties associated with appropriate track initialization.

Object Orientation Estimation

Further prior art relevant to the current concept focuses on the estimation of target orientation based on high-resolution radar, laser scanner or LiDAR measurements. These algorithms ignore the range rate measurement and focus exclusively on the locations of raw detections. FIG. 6 shows an example of the determination of target 9 orientation based on spatial distribution of raw detections. Example with three detections denoted with dots 6.

Several algorithms were proposed for such applications. They can be grouped in three main categories: Rotating Calipers based approaches; L-fit; and the Hough transform. Additional modifications to the rotating calipers method were proposed. The accuracies reported for these algorithms based on the assessment of experimental results were about 10 degree RMS. The accuracies were shown to be dependent on the relative target orientation in SCS (sometimes called aspect angle or exposition angle). A further simplified approach suitable for real-time executions has already been implemented for Short Range Radars. The enhancement proposed in that implementation is to treat the orientation of the tracked object associated with a considered cloud as an initial guess for instantaneous orientation estimation. Next the adjustment of the orientation is performed using the distribution of detections captured only in the current look. This approach can be treated as an optional ‘hot starting’ of the instantaneous orientation estimation. The results are promising in particular when the target is close to the host and a sufficient number of raw detections is available. It is noted here that the determination of target orientation in prior art is not related to the estimation of target motion. Two branches of building blocks were identified: development of algorithms considering only azimuth and range rate measurements and development of algorithms considering only azimuth and range measurements

DESCRIPTION OF THE PREFERRED EMBODIMENTS

The invention includes provides instantaneous estimation of target planar motion from detections of a single radar look of e.g. a high resolution Doppler radar with the possibility of detecting that the target has a yaw rate above a certain threshold. It is to be noted that the “raw data” from this single radar look provides the parameters of r_(i)—range (or radial distance), θ_(i)—azimuth angle, {dot over (r)}_(i)—raw range rate (or radial velocity) for each ith point of m point detections on a rigid body. These are the parameters which are used to determine the yaw rate as well as longitudinal and lateral velocities, where i=1, . . . , m. It is to be noted that the term instantaneous or single look radar data would include reflection data from a “chirp” in Doppler techniques which may scan over e.g. up to 2 ms. By this known methodology range rate may be determined.

Methodology according to aspects are based on the combination of the cloud algorithm and spatial (or geometry) based estimation of object heading. In this way all three radar measurements (i.e. range, azimuth and range rate) are exploited to overcome the previously known limitations.

For the sake of convenience, the earlier presented cloud algorithm for yawing targets is now expressed in sensor coordinate system. In an initial step, the method comprises emitting a radar signal at a single time-point instance and determining from a plurality (m) of point radar detections measurements therefrom captured from said radar sensor unit in a said single radar measurement instance the values for each point detection of range, azimuth and range rate; [r_(i), θ_(i), {dot over (r)}_(i)]. In the next step, estimates of range rate equation parameters (of the target) obtained by the Least Square Solution to:

${\overset{˙}{r}}_{i,{cmp}} = {{\left\lbrack {\cos\theta_{i}\ \sin\theta_{i}} \right\rbrack\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}.}$

Starting with the definition of over-the-ground velocity vectors at the location of each detection with longitudinal u_(t,i,scs) and lateral velocities v_(t,i,scs) aligned along sensor axes:

$\begin{bmatrix} u_{t,i,{scs}} \\ v_{t,i,{scs}} \end{bmatrix} = {\begin{bmatrix} {\left( {y_{t,{COR},{scs}} - y_{t,i,{scs}}} \right){\overset{\hat{}}{\omega}}_{t}} \\ {\left( {x_{t,i,{scs}} - x_{t,{COR},{scs}}} \right){\overset{\hat{}}{\omega}}_{t}} \end{bmatrix}.}$

The range rate equation becomes: {dot over (r)}_(i,cmp)=(y_(t,COR,scs)−y_(t,i,scs))ω_(t) cos θ_(i)+(x_(t,i,scs)−x_(t,COR,scs))ω_(t) sin θ_(i).

Since: y_(t,i,scs) cos θ_(i)=r_(t,i) sin θ_(i) cos θ_(i)=x_(t,i,scs) sin θ_(i) one can write: {dot over (r)}_(i,cmp)=(y_(t,COR,scs))ω_(t) cos θ_(i)+(−x_(t,COR,scs))ω_(t) sin θ_(i). The Least Square solution to this problem is: {tilde over (c)}_(t,scs)=y_(t,COR,scs)ω_(t), {tilde over (s)}_(t,scs)=−x_(t,COR,scs)ω_(t), and {tilde over (c)}_(t,scs)(y_(t,COR,scs))={tilde over (s)}_(t,scs)(−x_(t,COR,scs)). It is to be noted that the range measurement is not used in velocity estimation.

In the next step the orientation angle of the target (γ_(t,scs)) must be computed. This may be estimated from raw detection data such as (r_(i), θ_(i),) and may be estimated by L fit, Hough transform or rotating caliper methodology

In the next step the line l_(PH) perpendicular to the orientation of the target and passing through the center of the target (x_(t,c,scs), y_(t,c,scs)) is determined by:

${l_{PH}:y} = {{{{- \frac{1}{\tan\gamma_{t,{scs}}}}x} + {\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + y_{t,c,{scs}}} = {{a_{1}x} + {a_{0}.}}}$

In the next step, the line l_(CA) (equivalent to reference numeral 7 of FIG. 5 ) passing through the center of rotation and sensor mounting position is determined by applying the cloud algorithm:

${{l_{CA}:y} = {{{- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}x} = {{b_{1}x} + b_{0}}}},$

where b_(o) can be zero—so the term may be removed in another form of the equation.

In a fourth step the intersection point of the lines from the previous two steps are determined with:

${{\overset{\hat{}}{x}}_{t,{COR},{scs}} = \frac{b_{0} - a_{0}}{a_{1} - b_{1}}},{{{and}{\overset{\hat{}}{y}}_{t,{COR},{scs}}} = {\frac{{a_{1}b_{0}} - {a_{0}b_{1}}}{a_{1} - b_{1}}.}}$

It is to be noted that if both lines are parallel (i.e. α₁=b₁), then there is no solution and the yaw rate cannot be estimated. Additionally: if α₀≠b₀: the lines are parallel, there is no intersection thus no yaw rate, the target is moving along a straight line and where α₀=b₀: the lines are identical and nothing can be said about the motion of the target.

In step 5 if both the orientation estimation and cloud algorithm return plausible results, the yaw rate of the target calculate is calculated as:

${\hat{\omega}}_{t} = {{\frac{{\overset{\sim}{s}}_{t,{scs}}}{- {\overset{\hat{}}{x}}_{t,{COR},{scs}}}{or}{\overset{\hat{}}{\omega}}_{c}} = {\frac{{\overset{\sim}{c}}_{t,{scs}}}{{\overset{\hat{}}{y}}_{t,{COR},{scs}}}.}}$

It is to be noted that the estimated yaw rate depends on both the Least Square result for cloud algorithm as well as the position of COR. For example, for the same {tilde over (s)}_(t,scs) yaw rate will be higher when the center of rotation is calculated to be near the host vehicle and will be lower when the {circumflex over (x)}_(t,COR,scs) is further away from the host. FIG. 7 illustrates a geometrical interpretation of the method. Center of rotation estimated to be in the location of the dot 10.

Results

Results of estimating the yaw rate by the proposed technique was evaluated using experimental data. The experiments were collected for stationary host vehicle and a single target vehicle. Both vehicles were equipped with a differential GPS which accurately measures positions, velocities and yaw rates. In post-processing stage, all detections captured from the vicinity of the reference object were selected for calculating cloud algorithm. The orientation of the target was taken directly from reference data.

FIG. 8 shows an example of test results of a maneuver which comprised the target 9 driving on a circular track 12. FIG. 9 a to d illustrate the movement. The point of intersection of the two lines 13 and 14 (equivalent to lines lCA and lPH respectively) remains roughly at the same position, which is consistent with the actual maneuver performed. The only exception is when the two lines are almost parallel. This confirms analytical results.

FIG. 10 shows an example of test results of a maneuver where the target 9 is driving straight and FIG. 11 a to d illustrates temporal movement. The two lines 13 and 14 remain parallel. This is consistent with the expected behavior, because the target is not yawing.

FIG. 12 shows examples of test results of a maneuver where the target is maneuvering at a junction and FIG. 13 a to d illustrate temporal movement. The intersection point changes its position as the yaw rate of the target smoothly changes during the maneuver. This is consistent with expected behavior. Yet again, promising result for junction assist applications.

In the scenario discussed in this document, detections coming from a single sensor are considered. There are different ways detections from multiple sensors can be used when the target is in the common area of their FOVs: each sensor may estimate its own lines resulting from cloud algorithm and orientation estimation

In ideal case all lines should intersect at the same point, i.e. the center of rotation. Each sensor can calculate its own center of rotation. These points can then be used to determine resulting estimate for the center of rotation. Such processing should consider plausibility measures for sensor estimate. Detections from all sensors can be collected for improved orientation estimation, but each sensor can calculate a cloud algorithm result based on its own detections.

Using Filtered Heading

Estimation of object orientation usually require more detections than the execution of cloud algorithm. For the purpose of yaw rate estimation/detections, one may consider using target orientation from an associated tracker object. In such a case, yaw rate can be estimated based on only two raw detections. The yaw rate estimation is likely to be more robust, but it will not be an instantaneous result.

Include the Proposed Approach in Object Tracking Filter

If instantaneous estimation is not required, one can consider formulating an object tracking filter. The measurement model of such a filter would then include estimations of yaw rate and center of rotation as proposed in this work.

Propagation of Uncertainty

If the uncertainty of cloud algorithm, as well as uncertainty of orientation is known then uncertainty propagation rules can be used for uncertainty estimation of center of rotation and yaw rate. It can be the most efficient way to determinate if motion estimation is reliable. Estimated uncertainty can be also used in Kalman filter as a level of measurement noise. In its pure form, the approach described here is not based on tracking or filtering and therefore does not require any motion model and no initialization stage is required.

The approach has a lower number of restrictions in comparison to currently known concepts based on Cloud Algorithm. It does not need measurements from several time instances and works for stationary host vehicle. It does not need the target to be in the field of view of two sensors mounted significant distance away from each other. The approach does not require filtering/tracking and does not rely on any target motion model. The algorithm is suitable for application in production embedded systems because of its low computational complexity. The method can be immediately used in state-of-art short range radars for estimation of objects at low range (up to about 15 meters); e.g. it is suitable for junction assist applications (crossroads, roundabouts). The accuracy of the results can only be improved if more detections from moving objects were available; this is likely to happen as sensing technology progresses and in sensors dedicated to Autonomous Driving applications. 

What is claimed is:
 1. A method of determining, by a host vehicle equipped with a radar system, a yaw rate of a target vehicle in a horizontal plane, the radar system including a radar sensor unit configured to receive signals emitted from the host vehicle and reflected by the target vehicle, the method comprising: a) emitting a radar signal and determining, from a plurality of point radar detection measurements obtained by the radar sensor unit, range, azimuth, and range rate values for each point radar detection; b) determining, by using a Least Square method and in a sensor coordinate system, Least Square solution values of longitudinal and lateral components of a range rate equation of the target vehicle, where the longitudinal and lateral components of the range rate equation are based on the azimuth and range rate values, a longitudinal velocity of the host vehicle, and a lateral velocity of the host vehicle, an x-axis of the sensor coordinate system being perpendicular to the radar sensor unit and pointing away and ahead of the radar sensor unit, a y-axis of the sensor coordinate system being perpendicular to the x-axis and pointing right of the radar sensor unit; c) determining an orientation angle and a center of the target vehicle, the orientation angle being relative to the x-axis and a first line representing an orientation of the target vehicle, the center being based on the range and azimuth values and a second line that is perpendicular to the first line and passes through the center of the target vehicle; d) determining a position of a center of rotation of the target vehicle using a position of the radar sensor unit of the host vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; e) estimating the yaw rate based on the position of the center of rotation of the target vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; and f) operating, by a driving-assistance system of the host vehicle, the host vehicle based on the yaw rate of the target vehicle.
 2. The method of claim 1, wherein determining the position of the center of rotation of the target vehicle comprises: determining a third line that passes through the center of rotation of the target vehicle using the position of the radar sensor unit of the host vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; and responsive to a determination that the third line and the second line are not parallel, determining an intersection point of the third line and second line, the intersection point being the position of the center of rotation of the target vehicle.
 3. The method of claim 2, wherein the method further comprises: responsive to a determination that the third line and the second line are parallel, determining the yaw rate of the target vehicle cannot be estimated and that the target vehicle is moving along a straight line.
 4. The method of claim 1, wherein the method further comprises estimating the orientation angle using one of an L fit transform, a Hough transform, or a rotating caliper algorithm on the range and azimuth values.
 5. The method of claim 1, wherein the plurality of point radar detection measurements is obtained from a single frame of radar data.
 6. The method of claim 1, wherein the second line has a slope equal to a negative inverse of a tangent of the orientation angle of the target vehicle and a y-intercept based on the center of the target vehicle and the orientation angle of the target vehicle.
 7. The method of claim 6, wherein the third line has a slope equal to the negative quotient of the lateral component and the longitudinal component of the range rate equation and a y-intercept equal to zero.
 8. The method of claim 1, wherein the yaw rate is estimated as: the longitudinal component of the range rate equation divided by a y-coordinate of the position of the center of rotation of the target vehicle; or negative one times the lateral component of the range rate equation divided by an x-coordinate of the position of the center of rotation of the target vehicle.
 9. The method of claim 1, wherein the method further comprises: determining estimates of longitudinal velocity and lateral velocity of the target vehicle based on the yaw rate and coordinates of the center of rotation of the target.
 10. The method of claim 1, wherein the first line passes through the center of the target vehicle.
 11. A system of a host vehicle comprising: a radar system configured to determine a yaw rate of a target vehicle in a horizontal plane, the radar system configured to: a) emit a radar signal and determine, from a plurality of point radar detection measurements obtained by a radar sensor unit of the radar system, range, azimuth, and range rate values for each point radar detection, the radar signal being reflected by the target vehicle; b) determine, by using a Least Square method and in a sensor coordinate system, Least Square solution values of longitudinal and lateral components of a range rate equation of the target vehicle, where the longitudinal and lateral components of the range rate equation are based on the azimuth and range rate values, a longitudinal velocity of the host vehicle, and a lateral velocity of the host vehicle, an x-axis of the sensor coordinate system being perpendicular to the radar sensor unit and pointing away and ahead of the radar sensor unit, a y-axis of the sensor coordinate system being perpendicular to the x-axis and pointing right of the radar sensor unit; c) determine an orientation angle and a center of the target vehicle, the orientation angle being relative to the x-axis and a first line representing an orientation of the target vehicle, the center being based on the range and azimuth values and a second line that is perpendicular to the first line and passes through the center of the target vehicle; d) determine a position of a center of rotation of the target vehicle using a position of the radar sensor unit of the host vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; and e) estimate the yaw rate based on the position of the center of rotation of the target vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; and a driving-assistance system of the host vehicle configured to operate the host vehicle based on the yaw rate of the target vehicle.
 12. The system of claim 1, wherein the radar system is further configured to determine the position of the center of rotation of the target vehicle by: determining a third line that passes through the center of rotation of the target vehicle using the position of the radar sensor unit of the host vehicle and the Least Square solution values of the longitudinal and lateral components of the range rate equation; and responsive to a determination that the third line and the second line are not parallel, determining an intersection point of the third line and second line, the intersection point being the position of the center of rotation of the target vehicle.
 13. The system of claim 12, wherein the radar system is further configured to: responsive to a determination that the third line and the second line are parallel, determine the yaw rate of the target vehicle cannot be estimated and that the target vehicle is moving along a straight line.
 14. The system of claim 11, wherein the radar system is further configured to estimate the orientation angle using one of an L fit transform, a Hough transform, or a rotating caliper algorithm on the range and azimuth values.
 15. The system of claim 11, wherein the plurality of point radar detection measurements is obtained from a single frame of radar data.
 16. The system of claim 11, wherein the radar system is further configured to determine the second line as having a slope equal to a negative inverse of a tangent of the orientation angle of the target vehicle and a y-intercept being based on the center of the target vehicle and the orientation angle of the target vehicle.
 17. The system of claim 16, wherein the radar system is further configured to determine the third line as having a slope equal to the negative quotient of the lateral component and the longitudinal component of the range rate equation and a y-intercept equal to zero.
 18. The system of claim 11, wherein the radar system is further configured to estimate the yaw rate as: the longitudinal component of the range rate equation divided by a y-coordinate of the position of the center of rotation of the target vehicle; or negative one times the lateral component of the range rate equation divided by an x-coordinate of the position of the center of rotation of the target vehicle.
 19. The system of claim 11, wherein the radar system is further configured to: determine estimates of longitudinal velocity and lateral velocity of the target vehicle based on the yaw rate and coordinates of the center of rotation of the target.
 20. The system of claim 11, wherein: the radar system comprises multiple radar sensor units, each radar sensor unit collecting a respective plurality of point radar detection measurements; and the radar system is further configured to: determine, from each respective plurality of point radar detection measurements and for each radar sensor unit, respective range, azimuth, and range rate values; determine, for each radar sensor unit, respective Least Square solution values of longitudinal and lateral components of the range rate equation of the target vehicle; determine, for each radar sensor unit, respective orientation angles of the target vehicle, respective centers of the target vehicle, respective first lines, respective second lines, and respective third lines; responsive to determinations that the respective third lines and the respective second lines are not parallel, determine respective positions of the center of rotation of the target vehicle; and estimate the yaw rate based on the respective positions of the center of rotation of the target vehicle and the respective Least Square solution values of the longitudinal and lateral components of the range rate equation. 